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Abstract 

Kinetic Vlasov equation for collisional plasmas with BGK (Bhatnagar-Gross- 
Krook) collision integral is used. The case of arbitrary temperature (i.e. arbitrary 
degree of degeneration of electronic gas) is considered. 

From kinetic Vlasov equation the distribution function in square-law appro¬ 
ximation on size of transversal electromagnetic field is received. The formula 
for calculation electric current is deduced. This formula contains one-dimension 
quadrature. 

It is shown, that the nonlinearity account leads to occurrence the longitudinal 
electric current directed along the wave vector. This longitudinal current is perpen¬ 
dicular to the known transversal classical current, received at the linear analysis. 

When frequency of collisions tends to zero, all received results for collisional 
plasma pass in known corresponding formulas for collisionless plasma. 

The case of small values of wave number is considered. It is shown, that the 
received quantity of longitudinal current at tendency of frequency of collisions 
to zero also passes in known corresponding expression of current for collisionless 
plasmas. 

Graphic research of dimensionless density of the current depending on wave 
number, frequency of oscillations of electromagnetic field and frequency of electron 
collisions with plasma particles is carried out. 
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INTRODUCTION 


In the present work formnlas for calculation electric current in classical 
collisional Fermi—Dirac plasma are deduced. 

At the solution of the kinetic Vlasov equation we consider as in expansion 
distribution function, and in expansion of size of the self-consistent electro¬ 
magnetic held the quantities proportional to square of intensity of external 
electromagnetic held. 

At such nonlinear approach it appears, that an electric current pas two 
nonzero components. One component of electric current is directed along 
intensity of electric held. This electric current component is precisely the 
same, as well as in the linear analysis. It is "transversal" current. Thus, 
in linear approach we receive the known expression of transversal electric 
current. 

The second nonzero electric current component has the second order of 
inhnitesimality concerning intensity of size of electromagnetic held. The 
second electric current component is directed along the wave vector. This 
current is orthogonal to the hrst component. It is "longitudinal" current. 

Generating in plasma of the longitudinal current by the transversal 
electromagnetic held comes to light the nonlinear analysis of interactions 
of the electromagnetic held with plasma. 

Nonlinear ehects in plasma already long time H-ISI are studied. 

In works [T] and j3] nonlinear ehects in plasma are studied. In work [3] 
the nonlinear current was used, in particular, in questions of probability 
of decay processes. We will note, that in work [2] it is underlined existence 
of nonlinear current along a wave vector (see the formula (2.9) from j2]). 

Quantum plasma was studied in works 0- izzj. Collisional quantum 
plasma has started to be studied in work |[T3]. Then quantum collisional 
plasma was studied in our works [IS] - [IS] ■ In this works quantum collisional 
plasma with variable collision frequency was studied. In works [20] - [22] 
generating of longitudinal current by transversal electromagnetic held in 
classical and quantum Fermi—Dirac plasma pni, in Maxwellian plasma [21] 
and in degenerate plasma [22] was investigated. 
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In the present work formnlas for calcnlation electric cnrrent are dednced 
in classical collisionless plasma at any temperatnre (at the any degrees of 
degeneration of electronic gas). 


1 SOLUTION OF VLASOV EQUATION 


Let ns show, that in case of the classical plasma described by the Vlasov 
eqnation, the longitndinal cnrrent is generated and we will calcnlate its 
density. On existence of this cnrrent was specihed more half a centnry ago 
[2]. We take Vlasov eqnation describing behavionr of collisional plasmas 
with integral of collisions BGK (Bhatnagar, Gross and Krook) 

f + v| + e(EG[v.H])| = ,(;(0)_^), (1,1) 


In eqnation (1.1) / = /(r,v,t) is the distribntion electron fnnction of 
plasma, E, H are the components of electromagnetic held, c is the light 
velocity, p = mv is the electron momentnm, v is the electron velocity, u 
is the effective electron collision freqnency with plasma particles, = 
/eg(r, v) (eq = eqnilibrinm) is the local eqnilibrinm Fermi—Dirac distribn¬ 
tion. 


/eg(r,n) = 


1 + exp 


£ - ff(r ) 
ksT 


1 -1 


= [1 + exp(p2 - a(r))] ^ = /eg(r, P), 


£ = mv^/2 is the electron energy, /i(r) is the chemical potential of electron 


gas, kB is the Boltzmann constant, T is the plasma temperatnre, P = 
P/Pt is the dimensionless electron momentnm, px = mny, vt is the 
thermal electron velocity. 


Vx = 




is the dimensionless chemical potential. 


(kBT) 


kBT = = 


mv^ 
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is the thermal kinetic electron energy. 

Let ns consider, that in plasma there is an electromagnetic held, repre¬ 
senting the rnnning harmonions wave 


Let us consider, that the wave vector k is orthogonal to potential of the 
electromagnetic held, 

kA(r, t) = 0. 

Electric and magnetic helds are connected with the vector potential by 
equalities 

1 dA 

E =-—, H = rot A. 

c dt 

For dehniteness we will consider, that the wave vector is directed along 
an axis x, and the electric held is directed along an axis i.e. 


k = /c(l,0,0), E = £'„(x, t)(0,1,0). 


Hence, 


_ 1 dA iuj , 

= -ci« = 

C/c c/n 

H = —Ey ■ (0, 0,1), [v, H] = —Ey ■ (Uy, 0), 

UJ UJ 




kv. 


dpx 


+ (cj - /cn^) 


dpy\ 


and also 


xTi n u 

V, Hh:^ = 0, because ~ v. 

c/p c/p 


Let us consider linearization of locally equilibrium function of distribution 


where 


feq{P, x) = fo{P) + g{P)6a{x), 


/o{P) = [1 + 


a{x) = a + 5a{x)^ a = const, 

dh{P) _ 
da (1 -H 
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The equation (1.1) can be copied now in the following form 


at ox uj 


1 ( 1 , r 


= iyfo{P) ^ g{P)iy6a{x). (1.2) 

Size Sa{x) we will hnd from the law of preservation of particles number 

f(f f) - 0 

y Ue, n^27rhY 

From this conservation law we receive that 

'5«{n f 9(P)l^ = /[/ - /o{^)l- 


j (27rh)^ 

From this equation we obtain that 


{2-iThY 


5a{x) = 


[/ - /o(P)]d*P 

[ g{P)SP 


We notice that 

J g{P)(fP = 27r 

where 


dP 


OO 


dP 


1 + 


TT 


1 + 


= 7r/o(a), 


/o(a) = 


dP 


— (X) 


Therefore 


5a{x) = 


1 + eP^- 

-d— 

7rfo{a) . 


= 2 


dP 


l + e^^-«‘ 


[/-/o(P)]d3p. 


The equation (1.2) will be transformed now to the integrated equation 

-h (cj - kv^,)- -h 

Cc; L dp:, dpy\ 

1 


+ Vx^ + pf = pfo{P) - 

ot ox UJ 


+iyg{P) 


7Tfo{a 


[f - foiPWP 


(1.3) 
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Let us search for the solution of the equation (1.3) in the form 


/ = MP) + h + h, ( 1 . 4 ) 


where 

f2 ^ El ^ g2*(fcx-a.i)_ 

Let us operate with method consecutive approximations, considering as 
small parameter size of intensity of electric held. Then the equation (1.3) 
with help (1.4) equivalent to the following equations 


dfi , 9fi 

Vrr — -h lyjl 


dt 


dx 


eE,, 


UJ 


jdfo dfo 

-H (a; - kv^)^ 

'dpy 


+ ^9{P) 


7r/o(Q() 


hSP. (1.5) 


and 




UJ 


dt 

+ (^ - kvx) 


df2 ^/2 f 
+ Vx^ + l^j2 = 


dx 

dpy 


+ ^9{P) 


7Tfo{a 


f2d^P ( 1 . 6 ) 


From equation (1.5) we obtain that 

{ly -iuj + ikvx)fi = 


eE„ 


UJ 


kv, 


dfo 

' dp X 


+ (cj - kvx) 


dfo 


Here 


^1 = 


1 


7rfo{a) J 

Let us enter dimensionless parametres 


dpy 

fidfP 


+ ^9{P)M- 


91 = 


UJ 


TJ 


kxVT 


y 


kxVT 


k 

"= V 


(1.7) 
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jYiyrp 

Here q is the dimensionless wave number, kx = , is the thermal wave 

h 

number, Q is the dimensionless frequency of oscillations of the electromag¬ 
netic held. 

In the previous equation we will pass to dimensionless parametres 


eK, 


O.kxPT'^T 


i{qP^ - z)fi = 
QPx) 


Here 


We notice that 


Hence 


dP, 

z = Q + iy = 
dfo 


dP^ 


P. 




dPy 

(jj + iy 
kxVx 

dfo 
dP, 


+ y9{P)M 


( 1 . 8 ) 


-Py 


TD\ 

yPy-^ p - yPp- 


dP, 


= n 


dfo 

dP, 


Now from the equation (1.8) we hnd, that 


^ leEy dfo/dPy _ g{P) 

^ kxPrVT qPx- z qP:,- z 

Let us substitute (1.9) in the equation (1.7). We receive equality 


(1.9) 


1 + 


ly 


27r/o(a) 


5(P)^\ ^ leEy r dfo/dPy ^,^^ 

qPx- z j krPTVT J qPx - z 


It is easy to see, that integral in the right part of this equality is equal 
to zero. Hence 

Ai = 0 . 

Thus, it agree (1.9) function fi is constructed and dehned by equality 

leEy dfo/dPy 


fl = 


( 1 . 10 ) 


krPTVT qPx - z 

In the second approximation we substitute fi according (1.10) in the 
equation (1.6). 
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We will receive the equation 


le 




kVn 


V — 2 iuj + 2 ikvx)f 2 = 

d fdfo/dPy^ ^ ^ ^ ^ d / dfo/dPy \- 


krPrVT^ L qP^ - z 

+iyg{P)A2. 




dpy \ qPx - z /\ 


+ 


Here 


Ao — 


1 


hPP 


( 1 . 11 ) 


7r/o(a) 

Let us pass in this equation to dimensionless parameters. We receive 
the equation 


2i{qP^ -X- y )/2 = 


ie'^El 




+y9{P)A2. 


Let us designate 


, ^ iy (jj V (jjPiv/2 

y = n + — = --h L - ' 


2 kxVT 2kTVT 
From last equation it is found 


kxVT 


/2 = 


e^El 


2k^p^v^9l 


p d / dfo/dPy \ n - qPg, d^fo 

^ ^dPq,\qP^-z) qP,-zdP; 


1 


qPx - z' 


ly g{P) 


:Ao. 


( 1 . 12 ) 

2 qP,-z' ^ ^ ^ 

For hnding A 2 we will substitute (1.12) in (1.11). From the received 
relation it is found A 2 


Ao = 




Ji 


2ktp‘?pv‘^9l p / N , iy j 
^ TTfoia) + —Jo 


Here 


Jn = 


g{P)d^P 
qPx - z' ’ 
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T _ f \ u ^ f^fo/9Py\ , n-qP^d^fi 

' J ^dpSqP.-z)^ qP,-zdP^ 


d^P 


qPx - z' 

Substituting the found value A 2 in (1.12), definitively we find function 
/2 in the explicit form 


/2 = 


e^E?. 


2k‘^p^v‘^Q, 


d / dfo/dPy \ ^ n - qP:, d‘^f( 


dPx \qPx- zJ qPx- z dP^ 


1 


qPx - z' 


where 


+7 


7 


e^E?. 


g(P) 


2k‘^p^v‘^Q. qPx — z' ’ 

^ (^Z//2)Ji 

7r/o(a) + {iy/E)JQ 


(1.13) 


(1.14) 


2 DENSITY OF ELECTRIC CURRENT 


Let us hnd electric current density 


j 



2ddp 

{27Th)^ 


( 2 . 1 ) 


From equalities (1.4) - (1.6) it is visible, that the vector of current 
density has two nonzero components 


j {jx: jyi 0 ) • 


Here jy is the density of transversal current. 



2dz'p 


= e 


'^yh 


2d^p 

{2TihY 


This current is directed along an electric held, its density it is dehned 
only by the hrst approximation of distribution function. 

The second approximation of distribution function the contribution to 
current density does not bring. 

The density of transversal current is dehned by equality 



2ie^p\ 

{ 27 Th)^kT 



{df 0 /dPy)P, ^,^ 
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This current is proportional to the hrst degree of size of electric held 
intensity. 

For density of longitudinal current according to its dehnition it is had 

f f 2d^p 2evTPT /“ d r nS d 


(27rh)^ J (27rh)^ (27rh)^ 

By means of (1.6) from here it is received, that 


Jx 


e^Eyfri 

{27Thyk‘^Q 


d / dfo/dPy \ X - qP^ d‘^fo 

^ ^ dPx\qPx - z ) qPx-zdP^ 


Pxd^P 

qPx - z' 


+7 


PMP)d^P 


( 2 . 2 ) 


qPx - z' 

In integral from the second composed from square bracets (2.2) internal 
integral on Py it is equal to zero: 

Py = -\-00 


OO 


dpf >' ~ dP„ 


— OO 


= 0 . 


Py = — 00 


In the hrst integral from square bracets (2.2) internal integral on P^ is 
calculated in parts 


OO OO 

f d ( dfo/ap„ \ p^dp, j r 
J aPi \qP^- z) qPi- z' ^ J 


dfo/dPy 

(qPi - z' f'iqPx - z) 


— OO 


Hence, equality (2.2) becomes simpler 


Jx — 


e^EyTTi 

{27Th)^k^Q 


z q 


Py{dfo/aPy)<fP 
{qPy - z'YiqPy - z) 


+7 


Pyg[P)SP 
qPi - z’ 


Internal integral on variable Py we will integrate on parts 


Py^^<^Py = Pyfo 


— OO 


Py = + 00 OO 


Py— OO _0Q 


- j h(P)dPy = - j h(P)dPy. 


— OO 
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Thus, expression for the longitudinal current becomes 


Jx 


e^Eyin 

{27Thyk^Q 


+7 


z q 


MP)d^P 


(qPx - P f(qPx - x) 


PxgjP'jd^P 
qPx - P 


(2.3) 


Internal integral in plane (Py, Pz) we will calculate in the polar coordi¬ 
nates 


h(P)d?P 


oo 


dP^ 


OO OO 


{qP^ - z'Y{qP:^ - z) J {qPj: - z'YiqP^^ - z) 

ln(l -H e°‘~^^)dPx 

{qPx - z'YiqPx - z)' 


foiP)dPydPz = 


— OO 

00 


— 00 —00 


= TT 


— 00 


since 


00 00 



/o{P)rfP,rfP, = 7rln{l + e“-^0- 


— 00 —oo 


In addition 


qPx - P 


Pxg{P)S‘P f Pxh(Px)dPx 

qPx - z' J 

—oc 

rdr 


e“ rdr 


TT 


TT 


{1 + P'^~^){qT — z') J (1 + 

— OO 
OO 

Tiq f In(l -h )dT 


{qr — z'y 


— OO 


Equality (2.3) is reduced to one-dimensional integral 


Jx 


ire^EyTTiq 

{27ThYk‘yQ 


00 


00 


ln(l -|- e“ )dT 7 f ln(l + )dT 


[qr — z'Y{qT — z) 2 J [qr — z’^ 
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Let us rewrite the previous equality in the form 

.3 tp2, 


Jx 


TTe’^Eyinq 


z'Ju — ^'^02 


where 


Ji2 — 


{27TJi)^k‘^Q 

oo 

ln(l + )dT 


( 2 . 4 ) 


— CX) 


and 


Jo2 — 


{qPx - z){qT - z'Y' 


ln(l + e“ ^^)dT 
{qr — z'y 


— OO 


Let us return to consideration of size 7 . We will calculate integrals, 
entering in (1.14). We will calculate the hrst integral 


T = f Lp ^ ( 9fo/9Py \ , X - qPx d^fi 

' J ^ ^dP^qP,-z)^ qP,-zdPy^ 


d^P 


qPx - z' 


As it was already specihed, the integral from the second composed is 
equal to zero. The second integral as well as earlier we will calculate in 
parts. As a result n,e receive 

j =„ fp 5 P/oPJ’q <PP 

^ dP^\ qP^ - z J qP^ - z' 

^ 2 f Py[dfo/dP,]d^P 

^ J (qPi - ^)(qPi - py 

Now we will calculate in parts internal integral on the variable Py. As 
the result we receive 


Ji = -q^ 


h(P)SP 


(qP, - z)(qP, - z'Y' 
This integral has been calculated earlier. Hence 


Ji = - 7 Tq 


ln(l + e“ '^^)dT 
{qPx - z){qT - z'Y 




12 - 


—00 
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Let us calculate the second integral from (1.14). We have 


Jo = 


(X) CX) oo 

g{P)d^P f dp 


qPx - z' J qPx- z 

— OO —oo —00 


- j j g{P)dPydP,. 


The internal double integral is equal 


OO 00 

j J g(P)dPydP, = TT 


1 


1 + 


- = ^k{P)- 


— 00 —00 


Hence, the integral Jq is equal 


Jn = TT 


OO 


foir)dr 
qr — z' 


— 00 


So 


7^/o(«) + ^'^0 = 7T [ — —^/o(r)(ir. 

2 J qr - z' 


— 00 


Thus, the constant 7 is found 


7 


where 


00 


Jqi — 


W 2J12 
2^ Joi’ 

qr — ft 


qr — z 


y/o(r)dr. 


— 00 


Now (2.4) it is possible to present the formula in the form 

.3 tp2^ 


Jx 


TTe'^Eyfnq 

{27rhkk^Q 


o W?/ , 2^02 

''+2 + 4 '' 


Jl2- 


(2.4') 


Let us hnd numerical density (concentration) of plasma particles, cor¬ 
responding to Fermi—Dirac distribution 


A' = / UP) 


2(kp f P‘^dP j / \ 

{27Th)^ {27ThY J 1 + 27t‘^ ^ 

0 
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where kT is the thermal wave number, kr = 


mvT 


OO 


lo{a) = J ln(l + e“ )dT. 

0 

In expression before integral from (2.4) we will allocate the plasma 
(Langmuir) frequency 


' 4:7Te^N 


UJr. 


m 


and number density (concentration) and last we will express through 
thermal wave number. We will receive 


J 


^ kEl 


ykrPT^ 167r/o(«)f^ - 


iy iy 2 ^ 02 ] 


J 12 , 


(2.5) 


where 


Ldp hm. 

^ krVT mvl 


p 


is the dimensionless plasma frequency. 
Equality (2.5) we rewrite in the form 

■long _ 


Jx"° = J{kl,y,q)(TifrkE^, 


( 2 . 6 ) 


where ai^tr is the longitudinal-transversal conductivity, J{Q,y,q) is the 
dimensionless part of current. 


eUt 


eh / huj„ \ 2 


^Ltr — 


e f \ 

prkr Px inv^ J krPr ^ ^ 




1 


167r/o(Q;)n 
If to enter transversal held 


o , W , W 2'^02] 

"+2 + 4 " 


Jl2- 


Etr = E 


k(Ek) ^ ^ _ q(Eq) 


UJ 


kEtr = -1E,H 


k^ c 

then the equality (2.6) we will to write down in invariant form 

= J(n,y,q)ai,trkEi = J(n,y, q)a, 
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Remark. From the formula (2.5) (or from (2.6)) it is visible, that 
at ^ = 0 (or z/ = 0) i.e. when collisional plasma passes in collisionless 
(z —)■ Q, z' —> n), this formula in accuracy passes in the corresponding 
formula from our work [20] for collisionless plasmas 


■ long _ 
Jx 


CF\^trkE 


2 ^ 

^ 167r/o((n) 


CX) 



— CX) 


ln(l + e“ ^^)dT 
{qr — Uy{qr — Q) ’ 


Let us pass to consideration of the case of small values of wave number. 
From expression (2.5) at small values of wave number it is received 


■ long _ _ 7 p2 _ ^ 

Jx — criti./c£/ 

^ lbTTlo{a)ilzz 


OO 

- J ln(l + = 


SttQzz' Sttitiuj V uj 

Remark. At z/ = 0 from this formula in accuracy turns out the corres¬ 
ponding formula from [2D| for longitudinal current in the case of small 
values of wave number in collisionless plasma. 



3 CONCLUSION 

On Figs. 1 and 2 we will present behaviour of real (Fig. 1) and imaginary 
(Fig. 2) parts of density of the dimensionless longitudinal current at H = 
1 , y = 0.01 depending on dimensionless wave number q at various values of 
dimensionless chemical potential. At small and great values of parametre 
q curves 1,2 and 3 approach and become indiscernible. The real part 
has at hrst the minimum, and then the maximum. With growth of the 
dimensionless chemical potential the imaginary part of density of a current 
has one maximum. 

Further graphic research of size of density of the longitudinal current 
let us spend for the case of zero chemical potential: a = 0 (Figs. 3-6). 
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On Figs. 3 and 4 we will present behavionr of real (Fig. 3) and imaginary 
(Fig. 4) parts of density of the longitndinal cnrrent depending on dimensi¬ 
onless wave nnmbers q in the case = 1 at various values of dimensionless 
frequency of electron collisions. At small and at great values of dimension¬ 
less wave number curves 1,2 and 3 approach and become indiscernible. 

On Figs. 5 and 6 we will present behaviour of real (Fig. 5) and imaginary 
(Fig. 6) parts of density of the longitudinal current in dependence from 
dimensionless frequency of oscillations of the electromagnetic held Q in 
the case q = 0.3. At increase of dimensionless wave number q curves 1,2 
and 3 approach and at g' > 0.7 practically coincide. 

In the present work inhuence of nonlinear character of interactions of 
the electromagnetic held with the classical collisional plasma is considered. 

It has appeared, that presence of nonlinearity of the electromagnetic 
held leads to generating of the electric current, orthogonal to the held 
direction. 

Further authors purpose to consider a problem of the plasma oscillations 
and a problem about skin-effect with use square vector potential in expan¬ 
sion of distribution function. 
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Fig. 1. Real part of dimensionless density of longitudinal current, ^ = l,y = 0.01. 
Curves 1, 2, 3 correspond to values of dimensionless chemical potential a = —5, 0, +5. 



Fig. 2. Imaginary part of dimensionless density of longitudinal current, 
n = l,y = 0.01. Curves 1, 2, 3 correspond to values of dimensionless chemical 

potential a = —5, 0, +5. 
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Fig. 3. Real part of dimensionless density of longitudinal current, = 1, a = 0. 
Curves 1,2,3 correspond to values of dimensionless collision frequency 

y = 0.001,0.05,0.1. 



Fig. 4. Imaginary part of dimensionless density of longitudinal current, 17 = 1, a = 0. 
Curves 1,2,3 correspond to values of dimensionless collision frequency 

y = 0.001,0.05,0.1. 
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Fig. 5. Real part of dimensionless density of longitudinal current, y = 0.01, a = 0. 
Curves 1,2,3 correspond to values of dimensionless frequency of oscillations 
of electromagnetic field 11 = 0.3, 0.4, 0.5. 



Fig. 6. Imaginary part of dimensionless density of longitudinal current, a; = 1, a = 0. 
Curves 1,2,3 correspond to values of dimensionless frequency of oscillations of 
electromagnetic field 12 = 0.3, 0.4, 0.5. 
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